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Abstract 

The decoherence induced on a single qubit by its interaction with the environment is studied. 
The environment is modelled as a scalar two-level boson system that can go through either first 
order or continuous excited state quantum phase transitions, depending on the values of the control 
parameters. A mean field method based on the Tamm-Damkoff approximation is worked out in 
order to understand the observed behaviour of the decoherence. Only the continuous excited state 
phase transition produces a noticeable effect in the decoherence of the qubit. This is maximal when 
the system-environment coupling brings the environment to the critical point for the continuous 
phase transition. In this situation, the decoherence factor (or the fidelity) goes to zero with a finite 
size scaling power law. 
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I. INTRODUCTION 



Decoherence is the quantum phenomenon by which the coherence of a quantum system 
can be destroyed when it is put in contact with a large environment jl, 2|. The Schroedinger 
equation is a linear differential equation, consequently any linear combination of solutions 
is also a solution of the problem. Thus, a general possible quantum state is a superposition 
of quantum states. Nevertheless, such a state does not appear in the classical macroscopic 
world. The decoherence interpretation of quantum mechanics [l| claims that this is due to 
the interaction with the environment, which destroys the quantum correlations between the 
states of the system, making it to transite from a quantum superposition state to a classical- 
like mixture of states. Moreover, only a small set of states take part of the classical-like 
mixture; they are called pointer states 

The study of decoherence is important for several reasons: i) it might be responsible for 
the emergence of classical properties out of the underlying quantum nature of the physical 
systems, ii) it is a major problem for the construction of a quantum computer since it 
will produce the loss of the necessary quantum entanglement. Thus, both for fundamental 
reasons (i) and for practical purposes (ii) it is important to characterize the decoherence 
process and its effects on the physical properties of a quantum system. 

Along this line of study, it is important to address the issue of the effect produced in 
the coherence of a quantum state when the environment evolves between different quantum 
phases. There have been several works on the relation between decoherence and an environ- 
mental quantum phase transition [3 8]. Recently, we have presented a novel phenomenon 
in which the decoherence of the system suffers drammatic changes when the environment 
crosses an excited state quantum phase transition (ESQPT) {9]. An ESQPT is a nonanalytic 
evolution of the system as the control parameters in the Hamiltonian vary. It is similar to a 
ground state quantum phase transition but affecting to excited states. Correspondingly, an 
ESQPT can be classified in the thermodynamic limit as first order, when a crossing between 
two excited levels is present, or continuous, when the number of interacting levels is locally 
very large at an excited energy but without crossings. 

In Ref. [9 1 we presented briefly the case of a qubit in interaction with an environment 
modelled as a two-level boson system undergoing a continuous ESQPT. We used a particular 
simple Hamiltonian in terms of single control parameter to model the environment in order 
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to show the main effect. Here we present a more extensive study of a similar system including 
both first and second order ESQPT, and more general sets of parameters. Together with the 
exact evolution of the system, we present a simple mean field treatment. We show that the 
decoherence is maximal when the interaction of the system with the environment produces 
second order ESQPT, while no noticeable effects are observed in the case of a first order 
ESQPT. For the former case, a finite-size scaling analysis allows us to postulate that the 
fidelity goes to zero as soon as the interaction between system and environment is switched 
on. We also show that mean field treatment provides a good description for the decoherence 
of the small system, except around the critical points. 

The paper is structured as follows. In Sect. II, we present the model for the environment 
and study the phase diagram and its relation with the density of energy levels. We then 
discuss the interaction of the environment with the system. In Sect. Ill we show results 
for the decoherence factor. Both exact numerical results for large boson number, and an 
analytic mean field method with simple extensions of the Tamm-Dankoff approximation are 
presented. In Sect. IV, results for the decoherence factor in the case of continuous and first 
order ESQPT, including a finite size scaling study for the decoherence factor (or fidelity), 
are discussed. Finally, in Sect. V we summarize giving the main conclusions of this work. 

II. THE MODEL 

Following [3], we will consider our system composed by a spin 1/2 particle coupled to a 
spin environment by the Hamiltonian H$s, 

H ss = Is ®H s + |0> (0| <g> H\ + |1) (1| ® H Xl , (1) 

where |0) and |1) are the two components of the spin 1/2 system, and Ao, Ai the couplings 
of each component to the environment. The three terms Hg, H\ and H\ 1 act on the Hilbert 
space of the environment; therefore, it evolves with an effective Hamiltonian depending on 
the state of the central spin Hi = Hg + H\., i = 0,1. The term H\. makes it possible that 
the environment crosses a critical point as a consequence of the interaction with the central 
spin |3J]. 

Considering the initial state [^seify) = ( a |0) + b\l)) |£(0)), where |£(0)) is the initial 
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state of the environment, the evolved reduced density matrix of the system is 

p s (t) = Tr £ \V S e(t))(*se(t)\ (2) 
= |a| 2 |0)(0|+a6*r(t)|0)(l| 
+ a*br*(t) |1)(0| + |6| 2 |1}(1|. 

The off-diagonal terms of the density matrix are modulated by the decoherence factor 
r(t) which is the overlap between two states of the environment obtained by evolving the 
initial state |^(0)) with two different Hamiltonians, 

r (t) = (tf(0)| e iHot e- iHlt |tf(0)) . (3) 

If the environment is initially in the ground state of H , |0, g), the decoherence factor, 
up to an irrelevant phase factor, is 

r(t)^(0,g\e- iH ^\0,g). (4) 

This quantity has the same form as the Loschmidt echo or the fidelity, and it contains all 
the relevant information about the decoherence process. 

To be more specific, let us introduce as an environment a two-level boson system described 
by a generalized Lipkin Model, whose Hamiltonian is 

H £ = a n t - ^ (5) 
where the operators n t and Q w are defined as 

n t = tH, = sh + t j s + u th, (6) 

in terms of two species of scalar bosons s and t. a and u are two independent control 
parameters, and the total number of bosons iV = n s + n t is a conserved quantity. 

It is worth to mention that this two-level bosonic Hamiltonian is completely equivalent 
to an SU(2) spin Hamiltonian, with long-range spin exchange interaction. The equivalence 
is defined by the inverse Schwinger representation of the SU(2) generators 

S + = t*s = (S-)\ S z = ^(tH-Js), (7) 

where S represents the total spin of a chain of iV 1/2 spins. 

4 



A. Mean field theory for 



In order to study the phase diagram of the Hamiltonian (jSJ) as a function of the control 
parameters a and u, it is usual to rely on a coherent state of the form 

\N,0) = e^^y {st+/3 * t} |0) , (8) 

where |0) denotes the boson vacuum. The corresponding energy surface as a function of the 
variational parameter (3 is the expectation value of Hg (J5J) in the coherent state (jSJ) 



F(N 8) ( N ^\ H e\N,P) 



N- 



P 2 



-ha - 4 + A(3u(a - 1) + (3 2 [a + u 2 (a - 1)] }. (9) 



(1 + /3 2 ) 2 

Minimization of the energy with respect to (3, for given values of the control parameters 
a and u, gives the equilibrium value (3 e defining the phase of the system in the ground state. 
The value f3 e = corresponds to the symmetric phase, and (3 e ^ to the broken symmetry 
phase. 

This Hamiltonian has a second order Quantum Phase Transition (QPT) along the line 
U — 0, and a first order QPT for u ^0. In the latter, the critical point is defined as the 
situation in which the minimum in the symmetric phase and in the broken symmetry phase 
are degenerate and their energies are equal to zero. The study of the phase diagram has 



been done in several publications 10]. Here we summarize its main features. 



/3 = is always a stationary point. For u — 0, the solution with = is a maximum 
for a < 4/5, and becomes a minimum for a > 4/5. In the case of a = 4/5, (3 = is an 
inflection point, a = 4/5 is the point in which a minimum at [3 — starts to develop 
and defines the antispinodal line. 

For w / there exists a region where two minima, one spherical and one deformed, 
coexist. This region is defined by the point where the (3 = minimum appears 
(antispinodal point) and the point where the (3^0 minimum appears (spinodal point). 
The spinodal line is defined by the implicit equation, 

3a Af i-U^Y\, do) 



3a-4 B\ V A 




where A= (4-3 a + 2 (a- 1) u 2 ) 2 and B = 36tu 2 (a- 1) 2 . For example, for to = 1/V2, 
a ~ 0.822559. 

• In the coexistence region, the critical point is defined by the condition that both min- 
ima (spherical and deformed) are degenerate. At the critical point the two degenerated 
minima are at fl e = and e = u/2, and their energy is equal to zero. The critical 
line is therefore defined as 

4 + ^ 2 

ac= 5T^- 

For example, for to = l/y/2, a c = 9/11. 

• According to the previous analysis, for u ^ there appears a first-order phase tran- 
sition, while for oj = there is an isolated point of second-order phase transition at 
a = 4/5. In this case, antispinodal, spinodal and critical points collapse to a single 
point. 

In Fig. [T] we present a schematic view of the phase diagram for the environment Hamil- 
tonian (J5J) in the u — a plane. 

The Hamiltonian ([S]) also displays an Excited State Quantum Phase Transition (ESQPT), 
which is analogous to a standard quantum phase transition, but taking place at some excited 
critical energy E v of the system. We can distinguish between different kinds of ESQPTs. As 



it is stated in 



llj , in the thermodynamic limit a crossing of two levels at E — E c determines 




(a)w = (b)w = l/V2 



FIG. 2: Energy levels of the Hamiltonian ([5|) as a function of a for N = 50 and two different values 
of U). 



c 



a first order ESQPT, while if the number of interacting levels is locally large at E = E, 
but without real crossings, the ESQPT is continuous. As the entropy of a quantum system 
is related to its density of states, a relationship between an ESQPT and a standard phase 
iransition at a certain critical temperature can be established in the thermodynamic limit 
12]. These kinds of phase transitions have been identified in the Lipkin model fl3|. in 
the interacting boson model [14J], and in more general boson or fermion two-level pairing 
Hamiltonians (for a complete discussion, including a semiclassical analysis, see [15]). In all 
these cases, the ESQPT takes place beyond the critical value of the Hamiltonian control 
parameter, implying that the critical point moves from the ground state to an excited state. 

In Fig. [2] we show the energy eigenvalues of the environmental Hamiltonian (jSJ) with 
N = 50 bosons as a function of the control parameter a for u = in left panel, and 
to = 1/ y/2 in right panel. In both cases, we see for a < a c a collapse of several levels at 
E « 0. In the right panel (to = l/v2) we can also see a second critical curve for E < that 
divides the level diagram in two regions: one in which levels behave smoothly, and another 
in which the level density increases and some crossings are observed. 

One simple way to analyze the phase diagram is by means of the density of states. 
To obtain an analytical approximation for this quantity, one can start from a coherent 
state similar to (JE]), in which real parameter /3 is replaced by the complex parameter z = 
tan (0/2) exp (i£), in terms of which the energy is expressed as "H(0, ^) = (N, cf>, £| H \N, <fi, £). 
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FIG. 3: Density of states of the Hamiltonian (|5]) for N = 1000, a 
corresponds to uj = and the dotted line to w = 



0.5. The dashed line 



A good approximation for the density of states can be obtained by counting how many levels 
are there in an energy window dE, 

p(E) = ± j dZd<f>\J(<t>,0\6(n(<t>,0-E), (12) 

where | J{4>, £)| is the Jacobian of the transformation (0, £) — > (p, q), p and q are the canonical 
coordinates of the Hamiltonian, and M is a normalization constant. 

In Fig. [H]we show the density of levels of the environmental Hamiltonian (J5J), calculated 
by means of Eq. (|12|) . for N = 1000, a = 1/2 and the same values of u as in Fig. [2J As it 
can be seen, the collapse of levels at E — gives rise to a cusp singularity of p(E) for both 
u = and u = l/y/2. In the latter case, there also exists a jump in the density of states 
for a fixed value E < (E ~ —125 for this value of u) consistent with the energy spectra 
of Fig. [2j Although not shown, similar results are obtained for other values of a and u. In 
particular, the jump in the density of states at a certain value E < only appears for u > 0. 
Therefore, two different kinds of ESQPT exist in excited spectrum of Hamiltonian (jSJ). If 
we keep the terminology of thermodynamical phase transitions and we take the number of 
levels up to an energy E, N(E) = j dEp(E) as the analogue of the free energy F(N, T), we 
can conclude: (a) there exists a continuous A quantum phase transition at Ec — 0, for any 
value of parameter u; (b) there also exists a first-order quantum phase transition at some 
critical energy EP < if u > 0. 

To estimate the critical energies at which these quantum phase transitions take place, we 
can rely on the energy surface T-L((f), £). In Fig. H]we show T-L((j), £)/N in the thermodynamical 
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FIG. 4: (Color online). Energy surface 7i(4>,^)/N in the thermodynamical limit N — > oo for 
a = 1/2 and u = l/v2- Contour curves are drawn on the base of the figure (see text). 

limit N — > oo for a = 1/2 and to = l/y/2. The curves drawn in the base of the figure are 
contour curves for fixed values of the energy "H(0, £)/N = E. Gray curves (red online) 

(2) 

represent different values of E around E c for the continuous phase transition. The solid 
gray (red online) line represents the critical point E c = 0; this is the only value for which 
the contour curve is non-analytic. On the other hand, black curves (blue online) represent 
different values of E around the critical energy at which the first-order ESQPT takes place, 
E { c \ In this case, the critical value is the one at which the island around £ = 7r appears, 
that correspond to a local minimum in the energy surface. This entails the appearance of 
another region in the (0, £) plane for which the equation TL((f),^)/N = E has a solution, and 
consequently the density of states p(E) suddenly increases. 

B. Coupling to a single qubit 

Since we are interested in relating the phenomenon of decoherence in a single qubit with 
the structure of phases and critical regions in the environment as defined by the Hamiltonian 
(JTJ), we propose as a coupling Hamiltonian H\ i = Xin t . Choosing Aq = if the qubit is on 
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state |0) and Ai = A if the qubit is on the state |1), the effective environment Hamiltonian 
for each component of the system results into 



a 



a n t — 



N 



{a + X)n t 



QtQt 

1 — a 

N 



QtQt ■ 



(13) 
(14) 



This means that the qubit only interacts with the environment when it is on state |1). 

The system-environment coupling parameter A modifies the environment Hamiltonian. 
For certain values of a and A, this modification entails a crossing of the critical lines. Similar 
p_ were p r ev,o Usly anal^ed by S eve r a, a utllorS ffl, , udymg w hefc a q _ 
quench that drives the environment through a QPT implies some kind of universality in the 
decoherence process. 



Using the coherent state approach 10], it is straightforward to show that Hi goes through 
a ground state QPT at 



A* = (1 — a) (4 + u 2 ) - a (15) 

for a < a*. Therefore, if A > A* the quench makes the environment jump from one phase 
to the other. 

The main purpose of this paper to show that an ESQPT, instead of a ground state QPT, 
indeed produces dramatic consequences in the decoherence process. Using the coherent state 
approximation, it is straightforward to obtain that the coupling between the environment 
and the qubit entails an energy transfer in the former one, which is equal to 

AE(N, fi, A) = (N , fi\ \n t \N, 0) = ^Y^- (16) 

Therefore, the critical coupling A c which leads the environment to the critical energy E c is 

E(N,f3) + AE(N,(3,\ c ) = E c , (17) 

valid for both first order critical energy E^ and second order one E^\ In general, this is 
a trascendent equation, and therefore Ac 1 '' and A c 2 ^ (the A's corresponding to E^p and E^\ 
respectively) have to be obtained numerically. 
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III. CALCULATION OF THE DECOHERENCE FACTOR 



In order to calculate the decoherence factor (j3J) the expectation value of Hi ( |14p in the 
ground state |0, g) of H f[T5j) is needed. The decoupling of the complete system-environment 
Hamiltonian into the independent Hamiltonians Hq and H\ for each qubit state allows an 
exact diagonalization for large systems. In the following two subsections we will describe 
the exact formalism and make a comparison with mean field techniques supplemented with 
a Tamm-Dankoff approximation (TDA) treatment of the excited spectrum. 



A. Exact diagonalization 

A general Hamiltonian in terms of s and t bosons including up to two body terms is, 

H st = ath + b(t ] s + sH) + cthsH 

+ d{t ] st ] s + sHsH) + e{t ] sth + thsH) + fthth (18) 

where a, b, c, d, e and / are arbitrary parameters. 

Both Hamiltonians, H ffTHl) and Hi (TH1) . are particular cases of H st f iTHj) with the fol- 
lowing parameters, 

a — 1 
a = a + A — 2 — 

a — 1 

h = u ^r 

c = 2 a ~ 1 



d 



N 
a — 1 

N 





„ a — 1 


e = 








f = 


aCt — 1 






A = 


a — 1, 



(19) 

where A is an irrelevant global shift in energy. 

The exact diagonalization of the st Hamiltonian ffTSl) . and consequently of H and Hi, 

reduces to the diagonalization of a tridiagonal matrix in the basis 

t z t JV-i 

|M) = . 5 |0), (20) 
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where |0) is the boson vacuum and < I < N. Therefore, the dimension of the Hamiltonian 
matrix is d — N + 1. 

The relevant matrix elements are, 



(Nl\H st \Nl) = al + ft 2 + c/(l + N — I), 



(21) 
(22) 
(23) 



(Nl\H st \Nl + 1) = b^(N - I) (I + 1) + ey/(l + 1)(N - I) I 



(Nl\H st \Nl + 2) = dy/(l + 2)(Z + 1)^{N - l)(N -I -I), 



being all the others equal to zero. The diagonalization of the corresponding tridiagonal ma- 
trix can be done easily even for large N values, providing the exact results for the eigenen- 
ergies and eigenfunctions of Hq and H\ and consequently allowing to calculate numerically 



B. The Tamm-Dankoff approximation 

Before applying the exact diagonalization techniques to study the behavior of the decoher- 
ence as a fuction of the set of model parameters and particularly in relation to the quantum 
phase transitions (first and second order) in the ground (QPT) and excited states (ESQPT) 
of the environment, we will introduce an extension of the mean field approximation based 
on the TDA but including two phonon anharmonicities. 

Let us consider the condensate boson of the state jSJ) as a ground state deformed boson 
in a rotated basis. Since two Hamiltonians are involved, H and Hi, let us formulate the 
approximation for both in terms of a generic Hi (i = 0, 1). The variational parameter j3 in 
the condensate could be different for both Hamiltonians; therefore, the notation /3, (i = 0, 1) 
will be used to distinguish between both cases. With this notation, the deformed bosons (g 
and e) for Hi are related to the initial ones (s and t bosons) by 



r(t). 




(24) 



(25) 
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In terms of the deformed bosons the ground state and the first excited states are 

M> = ^(rI, 3 )"|o>, (26) 
l^)=yp=rl(rl„r- 1 |0). (27) 

In this framework, higher excited states can be constructed by directly replacing a ground 
state boson condensate by an excited (5 boson; this procedure is known as the Tamm-Dancoff 
approximation (TDA) method. In addition, with this basis is possible to write a diagonal 
Hamiltonian in terms of the new bosons. If only one body terms are included, 



Hi « (i,g\Hi\i,g)+ ( (i,e|^|i,e) - )rt e r iie (28) 

= s i ,o + A et rJ ie r iie , 

where E i>0 = (i,g\Hi\i,g) and A 6i = ((i,e\Hi\i,e) - (i, g\H\i, g)). 

The calculation for r(t) @ involves the Hq ground state and the Hamiltonian H%. Thus, 
it is necessary to relate the intrinsic bosons for H and Hi. The relation between both boson 
families is given by 

v 

V 

where this sum is for p = g and p = e, and the coefficients of the needed transformation are, 



^= ^W +ftft) - (31) 
f ^vmvm w °- 0l) - (32) 

With the preceding transformation it is possible to write |0, g) in terms of the H\ intrinsic 
bosons 

1/4 , \W 

\^9) = ^\f gg T{ 9 + f 9 eA,e) 10} ■ (33) 
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Using the binomial expansion of (|33|) is then straightforward to calculate the decoherence 
factor r(t) using the TDA basis up to an irrelevant phase factor 



N / \ 
k=0 ^ ' 



(34) 



where A ei = ((1, e|ifi|l, e) — (1, g\H%\l, g)). A more compact expression for r(t) can be 
obtained using the transformation, 



-iA ei tk -i(A sl /2)tN i(A ei /2)t(N-k) -t(A ei /2)ifc 



Therefore, the decoherence factor r(t) in the TDA reduces to 



(35) 



r(t) = e-^W {{f gg f e^^ + ifge) 2 e -*( A -i/ 2 )*) 
The matrix elements required for calculating r(t) are 



N 



(36) 



and 



(i,g\Hi\i, 



+ 



N{N — 1) 



(i,e\Hi\i,e) 



1 (1-/3 2 ) 2 
— (a + c + / - 26A) + c ), , X (iV - 1) 



1 + 



i + A 



+ 



4(iV -2^ + e (A-/?f)) 

A 2 ^ , nI A 



;i + a 2 ) 2 



+ (a + c + /)r ^(iV-l) +26^(^-1) 



+ 



(AT — 1)(N — 2) 



(37) 



((c + 2rf)/3 2 + 2e/3f + /A 4 ). (38) 



(1 + /3 2 ) 2 

A simple inspection reveals that decoherence factor r(t) fl36|) does not give a good ap- 
proximation of the exact results (see below and js]). The modulus of r(t) is 

N 



r(t)\ = (£) e^K/ 2 ) 4 + (/ 2 e ) e-K/ 2 )' " = \f ge f N 



.'/.-/ 



ye 



+ e 



A ei t 



(39) 



As a particular example, let us consider 0i = and /3o 7^ 0, that is, the situation in which 
the coupling of the qubit to the environment forces the environments to cross the phase 
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transition from the broken phase to the symmetric phase. In this situation 

f « = A (41) 

From these expressions, it is straightforward to obtain that \r(t)\ oscillates between 

|r(f)L« = 1; (42) 



K*)Lr 



/3o 2 -l 



/5 2 + l 



a? 



0, for N — > oo. (43) 



o 

Therefore, we can conclude that TDA approximation including just one phonon excita- 
tions does not account for the decay of the envelope of the decoherence factor reported in 
9] (see below for more details). This evidence suggests to go further within the spirit of 
TDA by including the anharmonicities of the two-phonon excitations. For this purpose it is 
needed to construct the states two TDA excitations as 



\i,e 2 ) = 4= / 1 (rjj 2 (rt )^~ 2 |Q> . (44) 
1 1 y/2y/(N -2)\ 

From this state we derive the diagonal part of the Hamiltonian as 



Hi « (i,g\Hi\i,g) + ((i,e\Hi\i,e) - (i,g\Hi\i,g)^Tl e Ti 

(i,e 2 \Hi\i,e 2 ) f . _ IrrL . v i {i,g\Hi\i,g)\ w t 



+ ' 2 '" ' - M^M + ™*«->»' j r;.,r:,|-,,r,, (45) 
= E i)0 + A ei r ie r i)e + fi ei r ie r ie r i!e r i)e , 

where E ifi = (i, g\Hi\i, g), A £i = ((i,e\Hi\i,e) - (i,g\Hi\i, g)), and Q ei = 
(Mmi^l _ (i, e \Hi\i,e) + MIhM>y 

With the preceding transformation ( 13*3*|) we can obtain the decoherence factor in the 
improved approximation up to an irrelevant phase factor 

r = E ( k ) (U m ' k) U 9 ef k e-^ k +^-W (46) 

In addition to (I37p and ( 138|) . the only needed matrix element for obtaining r(t) is Q ei 
which follows from (i, g\Hi\i, g) , (i, e\Hj\i, e), given above and 
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(i, e 2 \Hi\i, e 2 



+ 



+ 



+ 




{N {1 2 l%2 ^ (( C + ^ + + ■ 



+ 



(47) 



Inserting ( )37j) . ( 1381) and (I4T1) into ( 1461) we arrive to the final form of the decoherence factor 
r(t) within the extended TDA approximation. In this case, a semi-quantitative analysis as 
the previous one cannot be easily done. A comparison with exact numerical calculations is 
performed in next section. 

IV. RESULTS 

In this section we present the main features of evolution of the system described by ([T|) 
under the influence of the environment given by (E])- A brief report of the relationship 
between the decoherence in the qubit and the excited state quantum phase transitions in 
the environment was given in [9(. Here, we extend the analysis, comparing the numerical 
results with the Tamm-Dankoff approximation, and also facing the case of u 7^ 0, that was 
not considered in [9j. As two paradigmatic cases, we deal with the cases a = 1/2, and u = 
and cu = l/y/2. Different choices for the defining parameters of the model give rise to the 
same qualitative results. 

A. Decoherence factor for the Continuous ESQPT 

All the information about the decoherence process induced by the environment m 
the central qubit is encoded in the decoherence factor (|4j). As mentioned above, for the 
Hamiltonian we are using there is always a continuous ESQPT independently of the value 
of u>. In addition, for u^fl there also appears a first order ESQPT. In this subsection we 
will analyze the effect of the continuous ESQPT on the decoherence factor, while the effect 
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FIG. 5: (Color online) |r(i)| for a = 1/2 and five different values of A for two selections of the 
coupling system-environment parameter u, co = on the left and u = l/v2 on the right. In all 
cases N = 1000. Solid (black) lines correspond to the exact solution, and dashed (red) lines to the 
TDA calculation. 



of the first order ESQPT on the decoherence factor will be discussed in the next subsection. 

In Fig. [5] we show the results for the decoherence factor for N = 1000 bosons and 
a = 1/2, for two tu values, cu = (left panels) and to = l/y/2 (rigth panels). The objective 
of this figure is to show the effect of the continuous ESQPT on the decoherence factor. 
Solving Eq. (|T7|) for the value of a = 1/2, the continuous ESQPT (E^ = 0) takes place at 

(2) (2) / — 

Ac = 0.75 for u = (left panel) and at A c = 1.17 for cu = l/v2 (right panel). Several 
features deserve to be discussed. First of all, we can see that the TDA calculation works 
pretty well for small and large values of A. In particular, the shape of the envelope, which 

(2) 

remains unaffected by the increase of A for A ^> A c , is very well described by the TDA 
calculation (see panels for A = 4 and A = 8 in Fig. E]). Since this approximation mainly 
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relies on the position of the first and the second excited states of Hz, we can conclude 
that the information contained in the low energy spectrum is enough to have a good idea 
about the properties of the highest excited levels of the environmental Hamiltonian. Note 
that switching on the interaction between the central qubit and the environment entails an 
effective increase of the environmental energy roughly given by AE = (go\ H\(X) \qq) — Eq, 
and therefore a large value of A implies that the state of the environment jumps from the 
ground state to a mixed high-energy state. 

On the other hand, as it is clearly shown in the left panels corresponding to A = A c = 0.75 

(2) 

and A = A* = 1.5, and the right panels A = A c = 1.17 and A = A* = 1.75, the TDA 
calculations fails for intermediate values of A. These two values correspond to the critical 

(2) 

couplings A c and A*, corresponding to the excited state and the ground state quantum 
phase transitions, given by Eqs. flTTj) and ffT5l) respectively. The reason why the Tamm- 
Dankoff approximation does not work for these values is straightforward. The ESQPT 
entails a singularity in the energy spectrum far above the first excited state, which gives rise 
to the main contribution in the TDA calculation. On the other hand, the ground state QPT 
does not affect the decoherence suffered by the central qubit because the coupling A makes 
the environment to jump far above the critical point which entails a singularity in the gap 
between the ground and the first excited states. However, as the TDA calculation for r(t) 
strongly depends on this gap, it is spuriously affected by the QPT induced by the critical 
coupling A*. 

Finally, the best agreement between the Tamm-Dankoff approximation and the exact 

(2) 

calculation happens for A = 0.3, far below A c . Not only the envelope of the decoherence 
factor is well reproduced, but also the positions of the local maximum are well placed. In 
this case, the small coupling makes the environment to jump from the ground state to a 
mixed low-energy state. Therefore, it is reasonable to assume that the description provided 
by the TDA, which only takes in consideration the first excited state and a global measure 
of the anharmonicites of the spectrum, is a better approximation for small values of A. 

1. Analysis of the critical behavior of the decoherence at the continuous ESQPT 

As it is shown in Fig. [5J the decoherence of the central qubit behaves in a singular way 

(2) 

for a critical coupling A c , which makes the environment to jump to the critical energy 
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(a)w = (b)a = l/2 

FIG. 6: (Color online) 

r max i n function of the coupling A, for different values of a, ui, and N. In 
left panel u) = 0. Black lines represent the case a = 0; dark gray (red online) lines, a = 0.4; gray 
(magenta online), a = 0.6; and light gray (cyan online), a = 0.7. Solid lines represent N =10 000; 
dotted lines, N=2500; and dashed lines, N=600. In right panel, a = 1/2, and lines represent the 
cases uj = 0.2, oj = 1/2, uj = l/v2, and u = 1, with the same color code than right panel. Arrows 
show the critical coupling A c provided by Eq. ()17|) . 

Ec = 0. As the density of states in both cases u = and w^fl display the same critical 
behavior around this value (see Fig. [3]), also the same singular behavior for the decoherence 
is expected. 

In. Fig. [6] we show r max (X), defined as the second maximum of \r(t)\ (the first maximum 
is trivially \r(t = 0)| = 1). The left panel displays the case u = for several values of a, and 
the right panel the case a = 1/2, for several values of u> ^ (see caption for details). We can 
see that the behavior of this quantity is the same for u = and u^0. It evolves smoothly 

(2) 

and independently of the size of the system for values far from the critical coupling Ac , 
provided by Eq. ( ITTj) and shown in Tab. HI In a small region around A ~ Ac 2 '', r max becomes 
sharp, and the value of the minimum depends on the size of the system A^; the larger is the 

(2) 

system, the smaller is r max (\ c ). Therefore, for both to = and w / 0, the decoherence 
factor behaves in a critical way around A = A c where r max (\ c ) undergoes a dip towards 
zero which is sharper and deeper for larger values of N. 

We now investigate the thermodynamical limit, by performing a finite size scaling anal- 
ysis. The largest system that we could treat exactly has a size of around A^ = 10000; going 
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oj = 


a = 1/2 


a = a = 0.4 a = 0.6 a 


= 0.7 


OJ 


= 0.2 w 


= 0.5 u = l/y/2 oj = 1 


A c 2) = 2 Ai 2) = 1 A c 2) = 0.5 A c 2) 


= 0.25 


a c 2) 


= 0.83 A c 2) 


= 1.01 A c 2) = 1.17 A c 2) = 1.45 



TABLE I: Critical couplings Ac for the eight cases depicted in Fig. [6] 




FIG. 7: 

r max (Ac ) in function of the size of the environment N, in a double logarithmic scale. Left 
panel represents oj = 0; right panel a = 1/2. Squares represent the case a = (left) and oj = 0.2 
(right); circles, a = 0.4 (left) and oj = 0.5 (right); upper triangles, a = 0.6 (left) and oj = l/\/2 
(right); lower triangles, a = 0.7 (left) and oj = 1 (right). Straight lines represent the best fit to a 
power law r maa; (Ac ; ) = AN~~< . 



beyond this value is very difficult since for a complete calculation of r max (X c ) all the eigen- 
values and eigenvectors of the environmental Hamiltonian are needed. Starting with systems 
of N = 100, we analize the finite size scaling along two orders of magnitude. 

In Fig. 13 we show how r max (A c ) evolves with the size N of the environment, both for 
oj = and several values of a (left panel), and a = 1/2 and several values of oj ^ 0. In all 

(2) 

the cases, a power-law scaling r max (A c ) ~ iV~ 7 is observed, and therefore we can expect 
that r max (X^) — y in the thermodynamical limit iV — y oo. Nevertheless, subtle differences 
between verying a with oj = and varying oj with a = 1/2 are observed. The results for 
the exponent 7, shown in Tab. [TT1 are very close to the proposed 7 = 1/4 [9| for u; = 0. 
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u) = 





a = 




a = 0.4 




a = 0.6 




a = 0.7 


1 


= 0.247 ±0.003 


1 


= 0.248 ±0.003 


7 = 


0.248 ±0.001 


1 


= 0.245 ± 0.003 








a = 


1/2 










uj = 0.2 




w = 0.5 




UJ = l/y/2 




UJ = 1 


7 


= 0.255 ±0.006 


1 


= 0.259 ±0.003 


7 = 


0.264 ± 0.008 


7 


= 0.284 ±0.001 



TABLE II: Finite size scaling exponents 7 for the cases depicted in Fig. [7] 
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FIG. 8: \r(t)\ for a = 1/2, u = \/y/2 and three different values of A. In all cases N = 10000. 



However, the numerical estimates seem to increase for larger values of uj; in particular, for 
the case uj = 1, the result for exponent 7 is significatively larger than 7 = 1/4. 

B. Decoherence factor for the first order ESQPT 

For the case uj 7^ 0, the Hamiltonian considered produce, in addition to the continuous 
ESQPT studied in the preceding subsection, a first order ESQPT at energy EP . This 
critical energy can be estimated calculating the local minima in the energy surface T-L((j),^), 
as it is shown in Fig. HI Inserting this value in Eq. (JT7J) a critical coupling Ac is obtained. 
For the case a = 1/2 and uj = l/y/2 the first order EQSPT is obtained at « 1.05. 

In Fig. IB] we show the exact result for the decoherence factor \r(t)\ for a — 1/2, u — l/y/2, 
and three different values of A around A = A^ ~ 1.05. The most significative result is that 
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no trace of critical phenomena are observed in \r(t)\ - the shape of this magnitude is smooth 
around A = A c . Moreover, Fig. |U] confirms that r max (X) also behaves in a smooth an 
size-independent way. The conclusion is, thus, that the first-order ESQPT does not affect 
the decoherence induced in the central qubit. 

V. SUMMARY AND CONCLUSIONS 

The decoherence induced in a single qubit by its interaction with the environment, mod- 
elled as a scalar two-level boson model, is studied. The environment presents a quantum 
phase transition from symmetric to non-symetric phases at around a = 4/5, which can be 
first order (u 7^ 0) or second order (u> = 0). In the non-symmetric phase, the environment 
also presents excited state quantum phase transitions (ESQPTs): a second order one for 
any u value at = 0, and also a first order one for u 7^ at an energy Ec < 0. We have 
shown that the second order ESQPT affects dramatically the decoherence factor which goes 
rapidly to zero. A finite size scaling study shows that in that case the decoherence factor 
goes to zero at the critical point following a power law. On the other hand, the first order 
ESQPT does not affect the decoherence of the central qubit. 

We have also shown that a mean field treatment provides a good description of the 
decoherence factor r(t), except in the regions around the critical points. Therefore, more 
sophisticated approximations are needed to obtain an analytical description of the critical 
behavior of r(t), and, particulary, to estimate the critical exponent A. 
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